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Abstract. In the presence of spin-orbit coupling and inversion symmetry of 
the lateral confinement potential a single electron does not exhibit matrix Berry 
phases in quasi-two-dimensional semiconductor quantum dots. In such a system 
we investigate whether many-body correlation effects can lead to finite matrix 
Berry phases. We find that the transformation properties of many-electron 
wavefunctions under two-dimensional inversion operation do not allow finite 
matrix Berry phases. This effect is exact and is independent of the form 
of electron-electron interactions. On the other hand, quasi-two-dimensional 
semiconductor quantum dots with lateral confinement potential without inversion 
symmetry can have finite matrix Berry phases. We find that many-body quantum 
fluctuations can change matrix Berry phases significantly in such systems. 



1. Introduction 

Electron spins in two-dimensional semiconductors may be manipulated electrically [U 
EJ[3l[Il[5l[6l[7l[8l[9]. It is more challenging to control single or a few spins coherently 
in confined nano quantum dots [TOl [Til El 021 HH [HI US]. One way to perform 
such a coherent control electrically is based on matrix Berry phases[T71 [TH1 [T5] . 
There are several semiconductor nanosystems with spin-orbit coupling terms [201 HI] 
that exhibit matrix Berry phases: they include excitons[22]. CdSe nanocrystals[23j, 
acceptor states of p-type semiconductors [24] . and ring spin filters [25]. Recently it has 
been demonstrated theoretically that it is possible to control electrically electron spins 
of II- VI and III-V n-type semiconductor quantum dots [26] and rings [27] by exploiting 
matrix Berry phases. In these systems spin-orbit terms are invariant under time- 
reversal operation [28] and the discrete energy levels are doubly degenerate, and these 
properties are responsible for the generation of matrix Berry phases 29J. 

Coherent manipulation can be achieved by changing external parameters 
adiabatically in time. According to the theory of matrix Berry phase [18] the 
groundstate of a doubly degenerate Hilbert subspace of these II- VI and III-V n-type 
semiconductor quantum dots changes adiabatically in time as 



*(«) = ci(t)*(t) + c 2 (m{t), 



(i) 
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where and <E>(i) are the instantaneous degenerate single electron or many-body 
eigenstates (the overline in $(t) means time- reversal state of $(£)). The time 
dependent Schrodinger equation for the expansion coefficients C\ and C2 of Eq.(p} 
can be written as 

ihC v = ^ ' A VW C W v = 1, 2, (2) 

w 

where A vw = fcY^k(Ak)v,w^[t anc ^ ^ k are ^ ne adiabatic parameters labeled by k. 
The time evolution of the groundstate is governed by the 2x2 non-Abelian vector 
potentials (NAVPs) between the degenerate eigenstates 

d 

{A k ) v , w = i{$ v \ — (3) 

where $1 = $(£) and $ 2 = 

II- VI and III-V semiconductor quantum dots usually contain several electrons, 
and many-body effects may affect the matrix Berry phase. Formal expressions for 
many-body NAVPs can be derived including many body exchange and correlation 
effects [30 . Correlation effects are taken into account by writing many body eigenstates 
as a linear combination of single Slater determinant wavefunctions. When odd number 
of electrons are present the groundstates are doubly degenerate. At each time instant 
t they can be written as linear combinations of Slater determinant states 

M M 

1$) = m=f$=j2dim, (4) 

i i 

where T is time reversal operator and M is the number of instantaneous Slater 
determinant states in the linear combinations. These states are time reversal states 
of each other. (We have suppressed t in the quantities appearing in Eq.(f5]), and from 
now on we will do so unless explicitly written). The diagonal elements of the NAVPs 
are 

(^) 1 , 1 =i($|A|$) =i ^ c *^L + ^ c * c .( Bfc ). i) (5 ) 

* i,j 

where the elements of the NAVP between Slater determinant states are 

(flfcki=*<*i|^|*j>. (6) 

(The Slater determinant states in this expression can have different total confinement 
energies). It can be shown that if (B^ij is non-zero one can find single-electron 
wavefunctions 4> p and <f> q so that 

d 

(B k )i,j = (ofc) P ,g = i(^p\g^-m)- (7) 

When the single electron eigenstates (p p and <$> q belong to different energy shells (ak) p , q 
are called single electron inter-shell NAVPs [30], The other diagonal elements (Ak)2,2 
are given by Eq.© except that $ is replaced by $. The off-diagonal elements are 

(A*)^=i<*i^-ii>=<E<^ L +E c ^ ( Bk ^ ' ^ 

with (Afc) 2 .i = (Ak)i2- Within this approach one can use a Hartree-Fock 
approximation based on single Slater determinant groundstates, and show that fermion 
antisymmetry does not change the value of the matrix Berry phase. 
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These formal results have not been applied to investigate the interplay between 
matrix Berry phase and many-body correlations of II- VI and III-V semiconductor 
quantum dots. For example, not much is known about how the effects beyond Hartree- 
Fock approximation, i.e., correlation effects, change the matrix Berry phase. The total 
many-electron Hamiltonian consists of four terms: the kinetic part Hk , confinement 
potential Vc, spin-orbit terms H so , and electron-electron interactions V: 

H = H K + V C + H so + V. (9) 

Note that Vc may or may not be invariant under two-dimensional inversion operation. 
However, H so is not invariant under two-dimensional inversion operation, and, 
consequently, the total Hamiltonian H is not invariant under two-dimensional 
inversion operation, irrespective of the invariance of Vc- This implies that eigenstates 
of the total Hamiltonian H are not eigenstates of two-dimensional inversion operator. 
In the absence of many-body effects it can be shown that, when the lateral electric 
confinement potential has inversion symmetry, i.e., Vc is invariant under two- 
dimensional inversion operation, the matrix Berry phase is absent [26] . This is 
because off-diagonal elements of single electron intra-shell NAVPs, Eq. |3]) , may vanish. 
However, in the many-electron case the NAVPs may take finite values since they are 
related to the single electron inter-shell NAVPs, which can be non-zero, as can be 
seen from Eq.([7]). It is thus unclear whether the matrix Berry phase remains zero 
or not. In addition, it is not understood how many-body correlation effects change 
quantitatively the matrix Berry phase when a distortion potential breaks inversion 
symmetry of the lateral electric confinement potential. Such a quantitative estimate 
should be valuable in understanding experimental results of matrix Berry phases. 

In order to investigate these issues we use the formal results of Eqs.([4])-([8]). We 
have investigated the effect of many-body correlations, and have found that they do not 
induce a finite matrix Berry phase when the lateral confinement potential is invariant 
under two-dimensional inversion operation. This is an exact result. The main physics 
is that although there is coupling between different single electron energy levels, 
many-body correlation effects cancel this coupling. On the other hand, for lateral 
confinement potentials without inversion symmetry we find that many-body quantum 
fluctuations change the matrix Berry phase. In this case it is difficult to investigate 
exactly correlation effects. We have treated them in an approximation that includes 
a finite number M of many-body basis vectors, and have performed a numerical 
computation to estimate the effect of quantum fluctuations on the matrix Berry 
phase. Our approximate calculation shows that the effect of quantum fluctuations on 
the matrix Berry phase becomes more significant as the ratio between the Coulomb 
strength and the single electron energy spacing increases. The main results of our 
investigation may be tested experimentally in semiconductor dots, as we discuss in 
Sec.4. 

Our paper is organized as follows. In Sec. 2 we describe our model Hamiltonian. 
In Sec. 3 we compute the matrix Berry phase when the lateral confinement potential 
is not invariant under two-dimensional inversion operation. Discussions are given in 
Sec.4. 

2. Model Hamiltonian 

The total single electron Hamiltonian of a II- VI or III-V n-type semiconductor 
quantum dot contains an electric confinement potential and spin orbit coupling terms. 
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An electron with the effective mass m* of such a system can be described by the 
Hamiltonian [55] 

hs =h K + h R , 

h K = -^r + U(r) + V{z), 
Zm* 

hR = c R (a x k v - a y k x ) , (10) 
where the two-dimensional lateral confinement potential is 

U(f) = \m*ujlx 2 + \m*ioy + V p (x, y), (11) 

and the vertical confinement potential is V(z) (Here the two-dimensional coordinate 
is f = (x,y)). An electric field E is applied along the z-axis and electrons are confined 
in a triangular potential V(z), and it is assumed that only the lowest subband along 
the z-axis is occupied. Thus in our model quantum dots are effectively quasi-two 
dimensional. The Rashba constant changes when the electric field E is varied. 
The potential V p {r) — e'y perturbs the two dimensional harmonic potential with the 
strengths u> x and u y . This potential can be realized by applying a constant electric 
field along the y-axis and its strength e' is controlled by the magnitude of the applied 
electric field along the y-axis. The crucial point about U(r) and V(z) is that they can 
be changed electrically, which provides a means to control coherently electron spins. 
The Rashba spin orbit term[20 is /ir with Pauli spin matrices a x , y and a momentum 
operator k x = i (similarly with k y ). The Dresselhaus term can be also included, 
but since it does not change results qualitatively we omit it here. The Hamiltonian, 
Eq. (fT0|) , represents a simple model of the physical system, but it has all the correct 
symmetries. It is invariant under time reversal operation and each eigenenergy is 
doubly degenerate. The Hamiltonian is not invariant under two-dimensional inversion 
operation r — > — r since the Rashba spin-orbit term breaks inversion symmetry. In 
order to build up many-body wavefunctions we need to first construct single electron 
eigenstates. Each of these wavefunctions consists of the spin-up and down components: 

F l(f) \ - ( E mll c m „t(r>nf) \ 
Fi(r) ) ~ { c mln)l {r\m'n' |) ) ' {U) 

where |mn) are eigenstates of two-dimensional harmonic oscillators. 



3. Breaking of inversion symmetry, correlations, and matrix Berry phase 

When the lateral potential has two-dimensional inversion symmetry the effect of many- 
body correlations will not produce a finite value of the matrix Berry phase. This can 
be shown to be an exact result. (See Appendix B). It follows from the transformation 
properties of the wavefunctions under inversion operation. It should be stressed that 
although the lateral potential has inversion symmetry the total Hamiltonian does 
not. When the inversion symmetry of U(r) is broken many-body correlations will 
induce a finite value of the matrix Berry phase. It is not possible to compute this 
effect exactly unlike the case when U(f) has inversion symmetry. In a previous work a 
truncated single electron 4x4 Hamiltonian matrix was used [26 . However, many-body 
states built from these approximate single electron wavefunctions do not adequately 
describe many-body correlation effects. In this paper we find improved single electron 
wavefunctions, and use them to build many-body wavefunctions. Here we will 
compute the degenerate groundstates approximately by using a finite number of Slater 
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determinant basis states, i.e., using M — 4 in Eq.Q. This approximation should be 
valid as long as the single electron energy spacing is larger than or comparable to the 
characteristic Coulomb energy scale. 

3.1. Single electron Hamiltonian matrix 

We employ an improved approximation of a 6 x 6 truncated single electron Hamiltonian 
matrix, whose eigenvectors can be written, according to Eq. tfT!?]) . as 

|0) = ccuulOO T > + cd,i, t |0, 1, T > + co,a,t|0, 2, T > 
+ co,oJO,(U > + co,iJO,14 > + co,2,i|0. 2 .-i >■ 

(13) 

These expansion coefficients c mna of the basis states \mna) are eigenvectors of the 
6x6 Hamiltonian matrix 
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(14) 



= 3uj y . Then the 
Ei = 3tkjy, and 



The basis vectors of this Hamiltonian matrix are ordered as 1 00 f), 1 01 |), |02 |), 
etc. To reduce the number of independent external parameters we can set the ratio 
between the harmonic frequencies be a constant, for example, u) L 
energies of the two-dimensional harmonic oscillators arc Eq — Ihujyj 
Ei = Ahujy. The eigenvalues of this matrix are 

1 
3 
1 
6 
1 



Xo 



- cos | - 



Xo 



Eq 



Eq 



6 



vo sin 
y/3 sin 



(15) 



where 



El 



125, <5 
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E R , 



A = - 18E (Ep + El 



B = 



108 (E t 
A 



12<5) 3 - 11664^(5 2 , 



sin 9 



B 



(16) 

VA 2 + B 2 VA 2 + B 2 v ' 

Because of time reversal symmetry each eigenenergy is doubly degenerate. The doubly 
degenerate wavefunctions of the energy shells with eigenenergies X\, X 2 , X 3 are 
denoted by (<f)%, 02 ), (03: 04 ), and (05, 06), respectively. There is some arbitrariness in 
choosing these eigenstates since new eigenstates may be obtained by applying unitary 
transformation to the old set in each degenerate energy shell[TH]. We choose the 
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expansion coefficients of the first, second, and third pairs degenerate eigenstates of 
Eq.HU) as 

f 



and 



c(l) 
c(2) 

c(3) 
c(4) 

5(5) 
c(6) 



/3VT 
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f 

W 2 
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-*I,0,-l,a;,^,7f), 
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-<5 2 *,0,-l,a!,/3 2 *,7 2 *), 

"3, #3, 73, $3,0, 1) , 

-610,-1,(4,I%,t$). 



(17) 



(18) 



(19) 



Here c(p) denotes the expansion coefficients {c mna (p)} of p'th eigenstate. The 
quantities a p , /3 P , j p , 8 P are too complicated and lengthy to give here, however, they 
are all purely real or imaginary. This choice of the eigenstates simplifies the calculation 
of the many-body NAVPs. When U(r) has inversion symmetry, i.e., Ep = 0, then cm 
and 7i are zero. 



3.2. Many-body Hamiltonian matrix 

Using the previous results for single-electron wavcfunctions we construct many-body 
groundstates. Let us assume there are three electrons in the dot. In order to calculate 
the many-electron NAVPs we need to compute the expansion coefficients cj and di of 
Eq.Q. We include three single electron energy shells, each with double degeneracy. 
In our approximation we truncate the number of Slater determinant wavefunctions to 
four with the lowest total confinement energies. They are 

=a+a+a+|0), 
|*!) = a+a+a+\0) = |* 2 ), 

l* 3 ) = 44410), 

1*3) = 44410) = 1*4), (20) 

where af creates an electron in the i'th single electron eigenstate \(pi), given in 
Eqs. (ri7|) - (fT9|) . The vacuum state is |0). The truncated many-body Hamiltonian matrix 
is 



/ E A 
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-b 
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E B J 



(21) 



where the matrix elements are 

(* 3 \H\* 3 )=E B , 
{Vx\V\* 3 ) =a, 
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(*i|7|* 4 ) = b, 
(*i\V\V 2 ) =0 

(* 3 |F|* 4 )=0. (22) 

The matrix elements a and b contain Hartree and exchange contributions. All the 
quantities E A , Eb, a, b change when the adiabatic parameters change because they 
depend on the single electron wavefunctions, Eqs. (fT7]) - (fTl?l) . that are functions of the 
adiabatic parameters. The groundstate eigenenergy is doubly degenerate with the 
value 

Eg = I (Ea + Eb-D'A , (23) 



2 

where 

D=(E A -E B f + A{\a\ 2 + \b\ 2 ). (24) 
One of the doubly degenerate groundstates has the expansion coefficients 

(ci, C2, c 3 , c 4 ) = -j= {-bK, -a*K, 0, 1) , (25) 



where 

{-E A + E B + Dh) 

K= WFTW) ■ (26) 

The other degenerate state is obtained by taking time reversal of this state 

{d u d 2 ,d 3 ,d 4 ) = ^={aK,-b*K, -!,{)). (27) 
yivi 

Since c mna {p) of Eqs. (ri"7|) - (fT9)) are purely real or imaginary it follows from the 
expression for the Coulomb matrix elements, Eq. (|C.ip . that a is real and b is imaginary. 
Then from Eq. (f2T))) we also see that K is real. This implies that the expansion 
coefficients Ci and di are always purely real or imaginary. 

When inversion symmetry present in U(f) the results given in Eqs. (|25p and (|27[) 
simplify since 6 = 0, which follows from 

b= (^^H^i) = (13M16) - (13M16) 

+ (23M26) - (23|w|62) = 0, (28) 

where v is the Coulomb interactions between two electrons. (This is because, as 
can be seen from Eqs. (|A.4[) - (|A.6p . 03 ( r )* ''<t>&{ r ) , <M r )* ( t ) &{ r ) , and <Ps(r)*(^2(f) are odd 
functions of r) . In this case the expansion coefficients of the doubly degenerate many- 
body groundstates are given by 

(ci,c 2 ,c 3 ,c 4 ) = — ^= (0, -aK,0, 1) 



(29) 



and 



(d 1: d 2 ,d 3 ,d i ) = -^=(0*. 0,-1,0). 
v 



(30) 

Note that d^dz may not be zero, but {Bk)2,3, = {a>k)±,h — from Eq. (|A.8|) . 
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3.3. Correlation and matrix Berry phase 

We are now ready to calculate the matrix Berry phase. The strengths of the distortion 
potential and Rashba constant are 

E P = e'e y /V2 and E R = c R /V2£ y . (31) 

We choose these parameters as the adiabatic parameters: Ai = Ep and A2 = E R (The 
single electron Hamiltonian depends on them, see Eq. (TH|) . As explained in Sec. 2 these 
parameters can be controlled electrically). The adiabatic path is elliptic 

(\ 1 (t),\ 2 (t)) = (E R (t),E P (t)) = 

(E Rx + AE R cos(ujt),Ep iC + AE P sin(wi)). (32) 

We use the parameters E R>C = 0.5E , E PjC = Q.3E , AE R = 0.35£ , A£p = 0.2i£ , 
and to = Eq/10. The following steps are implemented consecutively in the computing 
the matrix Berry phase: 

(a) Single electron eigenvectors are evaluated numerically from Eqs.(fT7ll- (fT9|) . 

(b) Coulomb matrix elements are computed using Eqs. (|C.2|l . 

(c) Many electron eigenvectors are evaluated from Eqs. (|25|) and (|27p . 

(d) The many-electron NAVPs , given by Eqs.© and ([5]), are evaluated on the various 
point on the closed adiabatic path in the parameter space. For this purpose the 
expansion coefficients di are differentiated numerically. Also we use that the 
diagonal elements of the NAVPs between Slater determinants, Eq. ([6]) are zero 

(Bk)i,i = ^ = °> ( 33 ) 

pGocc. 

where the sum over p indicates a sum over single electron states that 
appear in the Slater determinant state This follows from the fact 

that c mn (p) is always purely real or imaginary, see Eqs. (fl7| - (fT9|) . Thus 

{a k )p,p = ^J^[Y, m n C *mn(p) c mn(p)} = for CaCh P' ( Not<3 that tlie narm0nic 

oscillator states \mn) do not depend on Afc). In addition it follows from the 
orthonormalization (*f f j|\I r j) = Sij and Eq.® that 

{Bk)i,j — (Bk)j,i for i 7^ j. (34) 
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Figure 1. f(t) as a function of t/T for k = 0.01,0.2,0.6,1. 
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Using these results we find that the NAVPs, Eqs.© and 

C 



A 1 




P 



P 




A 2 = 







), are off-diagonal: 

(35) 



with P, Q real. The time dependent Schrodinger equation for the expansion coefficients 
C\ and C 2 , given by Eq.@, can then be written as 



if{t)Ci 



where the function 

fit) 



P^ + Q dX2 



(36) 



(37) 



dt dt 

For each point (Ai(£), A 2 (<)) on the adiabatic path the off-diagonal elements of the 
NAVPs are evaluated numerically, as described in steps (a), (b), (c), and (d) above. 
As the strength of the Coulomb interaction k = -j- / Eq increases f{t) varies more 
significantly, see FigJT] The calculated matrix Berry phase for k < 1 is given as 
follows [3 lj 

Ci(T) \ f cos X isin X \ / d(0) 
C 2 (T) J \ ism X cos X J I C 2 (0) 



(38) 



where T is the period of the adiabatic cycle. The parameter x = J f(t)dt of Eq. 
is shown in Figj2] For k > 1 Slater determinant states with higher total confinement 
energies than those four we have used need to be included in the many-body basis 
set. This also implies that single electron states with higher energies than those six 
we have used must be included. 




0.0 0.2 0.4 0.6 0.8 1.0 
K 



Figure 2. The dependence of \ on re, where \ characterize the 2x2 matrix Berry 
phase and re is the ratio between the Coulomb energy scale and single electron 
levels spacing, re measures the strength of quantum fluctuations. 



As a check on the correctness of our numerical procedures we have verified 
numerically that the value of cos(x) is independent of the choice of the set of degenerate 
groundstates. (Although the elements of the matrix Berry phase depend on the the 
choice of the basis states the trace of it is independent of the basis states). By 
taking the limit of vanishing the strength of Coulomb interaction k — > in Eqs. pS)) 
and ([27]) we see that the degenerate groundstates are (01,02,03,04) = (a, (3,0,0) 
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and (di, d 2 , (^3, d^) = (—/?*, a*, 0, 0) (This limit is somewhat delicate since N\ and 
K diverge). We can also use another possible set for degenerate groundstates 
(ci, c 2 , c 3 , c 4 ) = (1,0,0,0) and (di, d 2l d 3 , d±) — (0,1,0,0). The elements of the 
NAVPs with respect to these new degenerate groundstates are (Ak)n = (-Bfc)n = 0, 
(^4/c)i2 = (Bk)i2 — (1/0)34- However, we find that the value of cos(x) is the same in 
these different sets of ground states. 

Using the computed matrix Berry phase, Eq. (|38j) . we now evaluate the single 
electron occupation numbers, which can be measured in tunneling experiments. 
Combining Eqs.fl]) and ^ we find that the many-body groundstate at each time 
instant is given by 

|*) = C 1 5^c 4 |* i )+C a ^d i |* i ). (39) 

i i 

The probabilities that a single-electron eigenstate p is occupied at t=0 and T are, 
respectively, 

f P (0) = l(Ci(0)c,(0) + C 2 (0K(0))|% 

i 

f P (T) = £ l(Ci (^)c,(0) + C 2 (T)d l (0))\ 2 ip , (40) 

i 

where f p = (^\at a p \^) , q(T) = Ci(0), and di(T) = di(0) (a and di are given in 
Eqs. (|2S|) and (|27|) ). If the single-electron eigenstate p is occupied (unoccupied) in the 
Slater determinant state \ we define 9i p = 1(0). At k = 1 we find for the occupation 
number of the third single electron level /3(0) = 0.0266 and fs{T) = 0.4494. For the 
fourth single electron level we find / 4 (0) = 0.9492 and / 4 (T) = 0.5264. The difference 
between / p (0) and f p (T) reflects the presence of a matrix Berry phase. It would 
be interesting to measure these differences in the single electron occupation numbers 
before and after an adiabatic cycle. 

4. Discussions 

The Hamiltonian of II- VI and III-V n-type semiconductor quantum dots with spin- 
orbit terms arc not invariant under two-dimensional inversion operator. Despite 
this, whether the lateral confinement potential is or is not invariant has important 
consequences on the matrix Berry phase. Our investigation shows that many-body 
correlation effects do not generate a matrix Berry phase when the confinement 
potential is invariant under two-dimensional parity operation. This is an exact result. 
It holds despite that the inter-shell single electron NAVPs couple different single 
electron energy levels. However, when the confinement potential is not invariant under 
parity operation our approximate calculation indicates that correlations can affect the 
matrix Berry phase significantly. 

Our results can be tested experimentally in self-assembled dots with wetting 
layers }3~2l [55] or in gated n-type semiconductor dots[31]. These quantum dots have 
several attractive features: The lateral shape of the dot can be distorted electrically 
to induce breaking of two-dimensional inversion symmetry. Moreover, the electron 
number can be varied from one to several electrons. These electric means for 
control provide excellent opportunities to test systematically the effect of many-body 
correlations. We have investigated quantitatively how quantum fluctuations affect the 
matrix Berry phase when the strength of Coulomb interaction is smaller or comparable 
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to the single electron level spacing. In self-assembled dots the characteristic scale of 
the single electron level spacing is 10 — 40meV, which is larger or comparable to the 
Coulomb energy scale of lOmeV. However, in gated semiconductor quantum dots[34 
the characteristic scale of the single electron level spacing is a few meV, which is smaller 
than the Coulomb energy scale. In order to obtain accurate results for these dots one 
needs to include a large number of Slater determinant basis states and single electron 
states. Nonetheless, even for these systems the matrix Berry phase should be absent 
when two-dimensional inversion symmetry is present, which should be experimentally 
testable. It should be noted that the matrix Berry phase depends on the geometric 
properties of an adiabatic path 35J. 
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(A.l) 



Appendix A. Single electron states and lateral inversion symmetry 

When the lateral confinement potential U (r) has inversion symmetry the single 
electron eigenstates simplify. For a given 2-fold degenerate energy shell we choose 
[3D] one of the eigenstates as 

where F e (f) and F (f) are even and odd functions of f. Although \<p) has even and odd 
spinor components it is not an eigenstate of the parity operator since the Hamiltonian 
is not invariant under two-dimensional inversion operation due to the Rashba term. 
(Note that any linear combination of the two states of Eq. (jA.lj) can also be chosen as 
a single electron basis state in the degenerate Hilbert subspace). 

We define a single- or many-electron wavcfunction to have a A-type property 
under parity operation if the spin-up part changes sign under parity operation: 

Fo(rl \ ( -FZ(r) 
F e (r) J { F* e {r) 

A wavcfunction has a B-type property under parity operation if the spin-down part 
changes sign under parity operation: 

' Fe(r) \ ( F* e {r) 
Fo^ J \ -F* {r) 

Each eigenstate \<j> p ) can be labeled by a subscript p. When p is odd the spin-up and 
-down components of the wavefunction are, respectively, odd and even functions of r. 
When p is even the odd and even properties are reversed. 

In order to include many-electron physics we need to fix single electron eigenstates 
of not only the first shell, but also of the second, third, and etc. energy shells. Here 
we choose them in the following specific order 

M <-» 



A - type. (A.2) 



B - type. (A.3) 
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I05> 

etc. 



,106 



(A.6) 



Note that the wavefunctions of a degenerate pair are chosen to be time-reversed states 
of each other. We have chosen the single electron wavefunctions <pi,<f>3,... to have A- 
type property, and (f>2, <f>4,, ... to have B-type property under parity operation, as shown 
in Fig lAll This particular choice simplifies the calculation of matrix Berry phases in 
the presence of many-body correlation effects. This corresponds to fixing a convenient 
'gauge', i.e., a single electron basis set. 



A-type 



3,c 



l.o 



F l,e 



B-type 



Figure Al. Each degenerate pair of eigenstates consists of A- and B-types. 
These two types of eigenstates are time-reversed states of each other. As the 
subscript i in <j>i increases the transformation properties of 4>i alternate between 
A- and B- types. 



The many-electron NAVPs contain single-electron NAVPs via Eqs.© and 0. 
So we need to understand first the properties of single-electron NAVPs. We can 
choose the adiabatic parameters as Ai = 2ftw x and A2 = ^ , where the lengths are 

£ x>1/ = yJh/m*uj Xy y. (The single electron Hamiltonian depends on them, see Eq. (|10|) ). 
The adiabatic constant Ai may be varied using the gate potential of the dot and A2 
may be varied by changing the electric field E along the z-axis. The single electron 
intra- shellNAVP elements [18] are i(4>p\^^\4> P ) j where <p p and <p p are degenerate single 
electron eigenstates. The NAVP elements between A and B or of B and A states can 
be shown to be zero [3D]. Since (j> p and 4> p are either of A and B or of B and A the 
intra-shell NAVP elements are zero. On the other hand, from the transformation 
properties of the eigenstates, given in Eqs. (|A.4j) - (|A.6|) . we can show that the single 
electron inter-shell NAVP elements are 

(ofc)p, g 7^ if p + q even, (A. 7) 

and 

{a k ) p , q = ifp + q odd, (A.8) 

where <j) p and <j> q belong to different energy shells. Note that different single electron 
eigenstates can be coupled through (ak) P . q if p + q is even. Thus the off-diagonal 
many-electron NAVPs, Eq.®, can be written in terms of non-zero inter-shell single 
electron NAVPs. Nonetheless it is possible to show that many-electron matrix Berry 
phase vanishes. 
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9i(ri)<p 3 (r 2 ) <P4(r 3 ) <Ps( r 4 > Ps( r 5 ) = 

Fl^P, ) F 3 , (r 2 ) F 4 , e (r, ) F 6 , e (r 4 ) F 8ie (r 5 ) v 
F^r, )F v (r 2 )F v (r 3 )F 6 , (r 4 )F SiO (i' 5 ) J 



even 
odd 



B-type 



Figure A2. A B-type Slater determinant is a sum of N\ terms. One of these 
terms is shown. 



Appendix B. Absence of matrix Berry phase and lateral inversion 
symmetry 

Many electron states can be written as a linear combination of Slater determinant 
states In the following we will choose l^i), ^3), ■■■ as A-type Slater determinant 
states, and |* 2 ) = r|*i),|* 4 ) = T|* 3 ),... as B-type Slater determinant states. 
The time-reversed state of A-type single electron wavefunctions are of B-type, and 
vice versa. The Slater determinant states are chosen in the order of increasing 
confinement energy (^^Hx +Vc\^i)- The total number of electrons N = Na + Nb is 
odd with the number of A-type single electron wavefunctions Na and of B-type Nb ■ 
As explained in Fig |A2l if Nb odd the Slater determinant transforms like a B-type i.e., 
spin-down part of the wavefunction changes sign. On the other hand, if Na odd |^) 
transforms like a A-type, i.e., spin-up part of the wavefunction changes sign. It can 
be shown that the NAVP between A and B Slater determinant states is zero 

(B^j =i(%(A)\-£-\* j (B)) =0. (B.l) 

This is because NAVPs between A and B single electron states are zero. 

We find that a correlated degenerate groundstate, |$) or |<&), is either A- or B- 
type. This is because the many-body Hamiltonian matrix element between A-type 
and B-type Slater determinant wavefunctions is zero, (^i(A)\H\^j(B)) = 0: If 1$) is 
A-type and |<f>) is B-type then 

|$) = d|*i(A))+C3|tf 3 (4)>+... 

|¥) = d 2 \<f 2 (B)) + <U\9 A {B)) + ... (B.2) 

We see from these results that, for a given index i, if an expansion coefficient Cj 
of one degenerate groundstate is zero then the expansion coefficient d t of the other 
time-reversed groundstate is non-zero, and vice versa. 

The off-diagonal elements of the many-body NAVPs, Eq.([5|), are zero. This can 
be shown as follows: according to Eq. (|B.2p . for each i, we have Cj = or di = 0, which 
implies that the first term of Eq.® is J2i c i§T^ = °- From Eq. ([B~2)) we see that 
when i j and c*dj is non-zero then ^ and are of A- and B-types, respectively. 
But this implies (Bk)i,j = 0, and the product c*dj(Bk)ij = 0. The second term of 
Eq.© is thus J2i j c idj(Bk)ij = 0. An explicit example of this is given below Eq. tfBT)]) . 
There is thus a delicate interplay between the many-body expansion coefficients c*dj 
and the elements of the NAVPs between Slater determinant states (Bk)ij. Since 
c i§JT = and c*dj(Bk)ij — the off-diagonal elements of the many body NAVPs 



Effect of many-body quantum fluctuations on matrix Berry phases 



14 



are zero: (^4/c)i.2 = 0. The matrix Berry phase is thus absent for doubly degenerate 
correlated states when inversion symmetry is present. This is true at any level of 
approximation represented by the number of Slater determinant states, M, included 
in Eq.(|3]). Therefore, this is an exact result valid for M — ► oo. 



Appendix C. Coulomb matrix elements 

The diagonal and off-diagonal matrix elements of H, Eq . (f2Tj) . depend on two-particle 
Coulomb matrix elements between single electron eigenstates p,q,r,s, that are given 
in Eqs.(Ul-C5), 



(pq\v\rs) = ^ P aja q a s c* npnp (p)c* mqn(i (q)c mrnr (r)c msns (s) 



m p , rriq , m r , m s , 



x (mpn p ,rri q n q \v\m r n r ,m s n s ) 

(C.l) 

where the Coulomb matrix elements between eigenstates of two-dimensional harmonic 
oscillator are 

(m p np,m q n q \v\m r n r ,m s n s ) — e 2 [ d 2 k — - (m p \e' lkxXl \m r ) (n p \e ikyV1 \n r ) 

J 2irk 

x (m q \e' ik ^\m s )(n q \e^ k ^\n s ) 

with 



(C.2) 



(me \m 



(^) 1/2 (%^) m - m 'e-^i™r m '(¥) K < m 



(^r) 1/2 (-^) m '- m e-^Z-'-™^) (m < m' 

(C 

and Lagurre polynomials (x). Similar expression can be found for (n\e tkyy \n'} with 
£ y replacing £ x . 
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